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ABSTRACT 



Aims. We investigate the change in the surface abundance of '^C during the evolution along the AGB, aiming to constrain third 
dredge-up models. 

Methods. High-resolution, near-infrared spectra of a sample of AGB stars in the LMC cluster NGC 1846 were obtained. A cluster 
sample ensures a high level of homogeneity with respect to age, metallicity, and distance. The C/O ratio and the ratio of '^C/''C were 
measured and compared with our evolutionary models. 

Results. For the first time, we show the evolution of the C/O and '^C/'^C ratios along a cluster AGB. Our findings allow us to check 
the reliability of the evolutionary models and, in particular, the efficiency of the third dredge up. The increase in both C/O and '^C/'^'C 
in the observed 0-rich stars is reproduced by the models well. However, the low carbon isotopic ratios of the two C-stars in our 
sample indicate the late occurrence of moderate extra mixing. The extra mixing affects the most luminous AGB stars and is capable 
of increasing the abundance of '^^C, while leaving unchanged the C/O ratio, which has been fixed by the cumulative action of several 
third dredge-up episodes. We find indications that the F abundance also increases along the AGB, supporting an in situ production of 
this element. 

Key words, stars: AGB and post- AGB - stars: variables - Globular clusters: NGC 1846 



1. Introduction 

To understand the role of low and intermediate mass stars within 
the cosmic matter cycle the efficiency of the third dredge up is 
a critical quantity. During the short but decisive asymptotic gi- 
ant branch (AGB) phase the third dredge up is responsible for 
mixing the burning products and the material produced via the 
s-process to the surface from where it can be ejected into the in- 
terstellar medium by stellar mass loss (see e.g. Busso et al. l 19991 
for a review). Current stellar evolution models (e.g. Straniero et 
al. [19971 Herwig[2000; Busso et al. |200il Stancliffe et al. |2004l 
Straniero et al. '2006; Karakas & Lattanzio r2007b provide quite 
detailed predictions for the surface abundance changes during 
the AGB phase. Different models agree qualitatively. The depen- 
dency of dredge up on various parameters like mass and metal- 
licity has been explored in the models (Iben & Renzini 119831 
Straniero et al. |2003l l. Quantitative checks of the model predic- 
tions for third dredge-up efficiency have been so far restricted 
to AGB stars in the solar neighborhood (e.g.Lebzelter & Hron 
12003 ' Busso et al.'2001) or to indirect methods like the study of 
the progeny of AGB stars or synthetic stellar evolution (e.g. van 
Eck et al. 2001 Marigo et al. 1999 ). However, for studying ongo- 
ing nucleosynthesis and the third dredge-up itself, observations 
of AGB stars are needed. 

Among the material dredged to the surface, '^C is of spe- 
cial interest as mixing of this element changes an oxygen-rich 
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Star (C/0<1) into a carbon-rich one (C/0>1). This leads to sig- 
nificant changes in the chemistry and, as a result, in the atmo- 
spheric structure and the mass loss properties of the star. The 
value of the C/O ratio is thus an indicator of the nucleosynthesis 
and mixing processes inside the star. For a few bright field AGB 
stars, measurements of the abundance of C and its isotopes exist 
(Lambert et al. |T986l Harris et al. [T987l Smith & Lambert |T990] l. 
Values for C/O range between 0.25 and 1.6 (see Fig. 9 of Smith 
& Lambert [T990] l. 

Comparing the findings from field stars with nucleosynthesis 
and mixing models is hampered by the rather large uncertainty in 
luminosity and mass, two very critical parameters for the mod- 
els. Clusters of stars offer an excellent possibility to investigate 
several questions of stellar astronomy. They provide samples of 
stars that are homogeneous in age and, in general, also in metal- 
licity located at the same distance. Thus it is possible to deter- 
mine their evolutionary status and their mass more accurately 
than for field stars. This provides an important advantage for 
comparison with models of stellar atmospheres and evolution. 

To investigate highly evolved stars on the AGB, globular 
clusters seem the best choice due to the number of potential tar- 
gets they include and due to the age range in which they are 
found. Globular clusters of the Milky Way are not good can- 
didates when studying the effect of the third dredge up. They 
are, indeed, too old so that their present generation of AGB 
stars should have a very low envelope mass, too low for the 
occurrence of a substantial dredge up. The Magellanic Clouds, 
however, contain a population of intermediate age clusters (e.g. 
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Girardi et al. l 19951 1 with AGB stars in ttie mass range 1.5 to 
2Mq. 

In this paper we present measurements of the C/O ratio and 
the isotopic ratio '^C/'^C for a sample of AGB stars in the cluster 
NGC 1846. In a previous paper (Lebzelter & Wood i2007i) . the 
variability of the AGB stars in this cluster was discussed. We 
refer to this paper for a recent overview of the literature on this 
cluster and will only repeat a short summary here, listing the 
values for the global parameters used in the present paper 

NGC 1846 is an intermediate age cluster belonging to the 
LMC. It has a metallicity of_[Fe/H]= -0.49 (Grocholski et 
al.|2P06 ). Lebzelter & Wood (12007b have determined the mass of 
the AGB stars to be close to 1.8 Mq resulting in a cluster age of 
1.4x10^ years. A number of AGB stars have been identified by 
Lloyd Evans (1980^, and we use his naming convention through- 
out this paper Near infrared measurements and estimates of the 
bolometric magnitudes of the AGB stars have been presented by 
Frogel et al. ( 1990 and references therein). No stars with a high 
mid-infrared excess, an indicator of a very high mass loss rate, 
have been found (Tanabe et al. 1998^ Lebzelter & Wood '2007'). 

The cluster velocity has been determined by various authors 
in the past. Schommer et al. ( |1992| l report a cluster velocity of 
240 km s"' derived from measurements of two individual stars 
(23 1 and 248 km s"' , respectively). This value coincides with the 
result by Freeman et al. ( 1983 ) and is in good agreement with the 
most recent study of individual stars in this cluster by Grocholski 
et al. (I2OO6I 1. who find a mean velocity of 235 km s ' from their 
17 cluster members. 
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Fig. 1. Observed spectral ranges. Various features are identified. 
A fit to the spectrum with reff=3550K, logg=0.25, C/O=0.30 
and a carbon isotopic ratio of 22 is given. In this fit, the F abun- 
dance has been reduced (see text). Due to the way the wave- 
length calibration was done (see text) the H band spectrum is at 
rest wavelength while the K band spectrum is shifted in wave- 
length according to the stellar velocity. 



2. Observations and data reduction 

High-resolution, near-infrared spectra of 12 AGB stars in NGC 
1846 from the list of Lloyd Evans (1980) were obtained with 
the Phoenix spectrograph (Hinkle et al. 1998) at Gemini South 
during six half-nights in December 2005. We observed all 12 
targets in the K band between 23620 and 23700 A; for 10 of 
them, spectra were also taken in the H band between 15540 and 
15587 A. Spectral resolution was set to approximately 50000. 
For most of the targets we achieved an S/N ratio of 65 or bet- 
ter We obained two to three observations per target and setting 
on different locations along the slit to allow for background sub- 
traction. Total integration times between 45 and 90 minutes per 
target and wavelength range were used. The spectra of each star 
were extracted and co-added using IRAF 

Wavelength calibration in the K band was done using telluric 
lines in the spectrum of an early type star that was also obtained. 
That spectrum was then used to correct the telluric features in 
the K band spectra of the program stars using the IRAF task 
telluric. The wavelength range in the H band was selected to be 
almost free of telluric lines to allow for a better comparison with 
synthetic spectra. Thus we used stellar OH lines present in the 
spectra of the 0-rich sample stars for the wavelength calibration 
in the H band. The C-star spectra were then calibrated with the 
solution derived from the 0-rich stars. 



3. Data analysis 

3.1. Contents of the observed wavelength ranges 

The observed spectral ranges in the H and K bands are shown 
in Fig. [1] for the star LE8, i.e. for a representative of the O-rich 
case. Several spectral features are identified. The H band range 
was selected to include a CO band head (3-0) and several lines 



of OH. A few atomic lines can also be seen. Finally a number of 
CN lines are present within the observed wavelength range. 

In the C-rich case, the //-band spectra are completely domi- 
nated by lines of CN and C2. The CO 3-0 band head is present, 
but its long wavelength side is strongly affected by neighboring 
features of the above-mentioned molecules. Lines of CO and CN 
are found in the observed wavelength range from the K band. 

3.2. Synthetic spectra 

The synthetic spectra are based on model atmospheres that were 
calculated with the COMARCS code, a modified version of the 
MARCS code (Gustafsson et al. 1975 J0rgensen et al. 1992)1. 
The temperature and, accordingly, the pressure stratification of 
the models are derived under the assumption of a spherical con- 
figuration in hydrostatic and local thermodynamic equilibrium 
(LTE). The LTE implies chemical equilibrium that is evaluated 
in both the model calculation procedure and spectral synthe- 
sis, in either case taking all major opacity contributors into ac- 
count such as H2O, TiO, CO, and CN in the oxygen-rich case or 
HCN, C2H2, C2, etc. in the carbon-rich case. For stars showing 
only small amplitude variability, LTE turns out to be a valid ap- 
proximation, while potential uncertainties in the determination 
of abundances are due to a deviation from hydrostatic equilib- 
rium rather than to non-LTE effects. COMARCS makes use of 
updated atomic and molecular opacities calculated with COMA 
(Aringer 2000). In calculating the model spectra, spherical ra- 
diative transfer routines originating from Windsteig et al. ( 199^ 
were used. The atomic line data were taken from VALD (Kupka 
et al. I2OOOI 1. whereas the molecular data were compiled from 
several sources (see Cristallo et al. 120071 for an overview and 
Lederer & Aringer |2007| for more recently added line lists). 
Model atmospheres and synthetic spectra that were calculated 
following the method described above have been shown to de- 
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scribe the spectra of cool giant stars appropriately (e. g. Loidl et 
al. |200n Aringer et al. l20U2] l. 

We want to point out that we calculated the model structure 
and the synthetic spectra in a fully consistent way with respect 
to spherical geometry and the atomic and molecular input data. 
Each atmospheric model is determined by six parameters: mass, 
metallicity, effective temperature, surface gravity, C/O ratio, and 
microturbulent velocity. From the full set of parameters we kept 
the mass (M - 1.5 Mq) and the metallicity ([M/H]=-0.4) fixed 
(cf. next section for more details and references). The effec- 
tive temperature was varied in steps of 50 K, and the values for 
Tgff range from 2600 K to 3850 K in our model set. In terms 
of the surface gravity, we chose a stepsize of 0.25 in logarith- 
mic scale, i.e. log (^[cto i"^]), ranging from -0.50 to -1-0.50. 
The microturbulent velocity in the model calculations was set 
to ^ = 2.5 km s ' apart from the model for the star LEI 3 where 
we adopted a value of ^ = 3.5 km s^^ (see below for details). 
The solar abundances were taken from Grevesse & Noels ( |1993l l 
whereby an inherent value of C/O = 0.48 emerges. As for the 
variation in the C/O ratio, we altered the abundance of carbon 
to arrive at the desired value and left the other abundances un- 
touched. We changed the C/O ratio in steps of 0. 1 in general, but 
also reduced the stepsize to 0.05 when necessary for the fitting 
procedure (see next section). 

The wavelength range covered in the calculation of the syn- 
thetic specti-a was 6400-6450 cm"' (15504-15625 A) in the H 
band and 4215-4270 cm"' (23420-23725 A) in the K band. The 
synthetic spectra were calculated with a resolution of 300000 
and then convolved with a Gaussian to get a resolution of 
50000, which is about the value attained in our observed spectra. 
Subsequently, the macroturbulent velocity was accounted for by 
another convolution with a Gaussian. In the course of our anal- 
ysis we also varied the carbon isotopic ratio '^C/'^C. The solar 
value as given in Anders & Grevesse j 19891 ) is about 89.9. The 
variation in the isotopic ratio has de facto no effect on the struc- 
ture of the model atmosphere and thus was only considered in 
the spectral synthesis calculations. In the next section we give 
details in the variation of the above-mentioned parameters for 
the individual tai^gets. 

3.3. Abundance determination 

To determine elemental abundances we used a two-step pro- 
cess. As a starting point we used T^ff and L values derived from 
broad-band, near-infrared photometry as described in Lebzelter 
& Wood (2007). Assuming a stellar mass of 1.5 M0, we derived 
a starting value for the stellar radius and the surface gravity. 
Around these parameters we calculated a small grid of synthetic 
H band spectra and selected the best fit with the observed spec- 
trum both by visual inspection and by a approach. Several 
relative line strengths could be used to estimate the correctness 
of the fit. We varied T^ff, log g, and C/O ratio; among these three 
parameters, changes in log g had the weakest effect. Thus a small 
uncertainty in the mass or radius estimate of the star only has a 
minor influence on the derived abundance^. In the O-rich case, 
changes in Teff have the strongest effect on the OH/Fe blend at 
15570 A, but an increase in temperature leads to a general in- 
crease in the strength of all features. An increase in C/O makes 
the CO band head stronger and all the OH lines weaker at the 
same time. Additionally such a change leads to more prominent 

' We thus do also not have to bother about the small difference in 
mass between the one derived in Lebzelter & Wood ( I2007t and the value 
assumed for calculating the synthetic spectra. 



lines of CN. In the chosen wavelength range, the model can be 
fitted to a well-defined pseudocontinuum. All these facts allow 
us to make good fits of the spectra in the O-rich case as illus- 
trated in Fig.[T]and to derive a C/O ratio from that. To achieve a 
good match of the lines, the macroturbulence was set to 3 km s"' . 
The influence of other values was tested. For the fit of LEI 6 
and LEO, a macroturbulent velocity of 4 and 7kms"', respec- 
tively, was necessary to account for the broader lines observed in 
these stars. For LEO we also increased the microturbulence to 
3.5 km s"' as mentioned before. We note that it was not our aim 
to measure the metallicity of the stars from our spectra because 
the selected wavelength ranges include only a small number of 
weak or strongly blended metal lines. The chosen metallicity of 
[Fe/H]= -0.4 for our models is close to the literature value and 
the metallic lines in the synthetic spectra are in reasonable agree- 
ment with the observed spectra. The fitting of the C-stars turned 
out to be much more complicated (see below). 

In a second step, the results from the H band spectrum were 
used as starting parameters for fitting the K band spectrum. The 
latter is - especially for the O-rich case - rather insensitive to 
changes in Teff , log g, and C/O ratio. Thus these parameters can- 
not be determined from the K band spectra alone. The chosen 
value of the macroturbulence could be confirmed. The main pur- 
pose for obtaining and analyzing the K band spectra was to de- 
termine the '^C/'^C isotopic ratio, which with the C/O ratio is 
also a very reliable indicator of third dredge up. For this the O- 
rich spectra offer two (almost) unblended lines of '^CO (2-0 R18 
and 2-0 R19) and two blends including contributions from '^CO 
lines. The best fit was selected by minimizing the squared differ- 
ence in the spectral range around the two unblended '^CO lines. 
Again, fitting of the C-stars was more problematic and is briefly 
discussed separately below. 

Uncertainties of the abundances were estimated from the 
range of parameters of the model spectra that still gave a good 
fit. For the C/O ratio, the Teff , and the log g, value this was done 
from the fit of the //-band spectrum, for '^C/'^C from the fit of 
the /T-band spectrum. Because the carbon isotopic ratio has a 
negligible effect on the fit of the //-band spectrum, as verified 
by test calculations, the fitting procedure was not repeated with 
the determined '^C/'^C ratio. 

The K band spectra also include a blend of the '^CO 2-0 R21 
line with an HE line. When comparing this line blend with the 
synthetic spectrum using the isotopic ratio derived from the two 
unblended '^CO lines, it turned out that this blend could not be 
fitted properly. An improvement in the fit can be either achieved 
by a change in T^ff or by a modification of the F abundance. 
The first one was only possible in a very limited range due to 
the constraints set by the fit of the H band spectrum. Thus we 
modified the abundance of F until we achieved a good fit. This 
was only done for the O-rich stars because in the C-stars the 
blend is additionally affected by a nearby CN line. 

The lowest excitation lines of CO in our K band spectra 
could not be fitted well, especially in the case of the more lu- 
minous O-rich stars in our sample. An explanation for this could 
not be found; however, it may indicate a problem of the model 
structure in the outermost layers of the star, which may reflect 
dynamic effects or an outflow of material expected for variable 
stars. 

3.4. Carbon stars 

When fitting the C-rich stars, we were confronted with two ma- 
jor problems. The first one concerns the uncertainties of the line 
lists for the C-bearing molecules. The line lists of CN and C2, 
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both critical for fitting the //-band spectrum, produced a num- 
ber of features in the synthetic spectra that were obviously either 
at the wrong wavelength or of the wrong intensity. The CN list 
could be improved by replacing the wavelengths in the line list 
with the observed ones taken from Davis et al. (2005). This was 
possible for a number of strong lines, while the weak lines typi- 
cally had no counterpart in the observed catalog. After that, a few 
CN lines with wrong line strengths still showed up in the syn- 
thetic spectrum. One of these lines, located at 6417.405 cm"', 
was removed by hand from the CN list for the later calculation. 
This line is also not in the catalog of observed lines. 

No observed reference is available for the C2 line list. Thus 
correcting this list is only possible by visual comparison of our 
spectra of C-stars in NGC 1846 with synthetic spectra. Within 
the studied range in the //-band, 5 lines were removed com- 
pletely and 9 were shifted in wavelength to fit the observed spec- 
tra better. In this wavelength range, lines of the two molecules 
including a '■'C atom (i.e. '■'C'^C and '■'C'^^N) play an impor- 
tant role for some of the blends. It turned out to be impossible to 
check the quality of these line data or to improve the line lists in 
that case because determining the carbon isotopic ratio was one 
of the aims of our project. We come back to this point in 14. 31 

The second major problem was the absence of reliable 'con- 
tinuum' points in the C-star spectra. The main effect of decreas- 
ing Teff or increasing the C/O ratio is the amount of depression of 
the (pseudo)continuum by a large number of weak and densely 
spaced lines over the whole spectral range. Thus the overall ob- 
served line depths can be fitted with a variety of combinations of 
these two parameters simply using different scaling factors be- 
tween the observed and the calculated spectrum. At first glance, 
a possible point for the pseudocontinuum seems to be present 
close to 15571 A. However, closer inspection reveals an absorp- 
tion feature at that position in the observed spectrum that is miss- 
ing in the synthetic spectrum. Therefore we think that this point 
is also not useable as a reference point. Thus we are limited to 
look for the comparably weak effect of using the ratios of fea- 
tures that seem to rely on one or the other of the two parameters. 
This, however, is hampered both by the uncertain line lists as 
discussed above and by the fact that a deeper CN or C2 line can 
either be achieved by an increase in C/O or a decrease in tem- 
perature. For determining the C/O ratio, we thus focused on the 
CO 3-0 band head that seems to be the only feature decoupled 
from this trend, but a larger uncertainty of the result compared 
to the O-rich case was unavoidable. The CO band head weakens 
relative to its neighboring features as C/O is increased, and it is 
almost insensitive to temperature changes. 

Fitting the ^T-band is eased by the lack of strong C2 lines, 
but as in the O-rich case, the depth of the various spectral fea- 
tures is rather insensitive to changes in the global parameters. 
Unfortunately, one of the two unblended '^CO lines from the 
O-rich case is heavily blended by CN in the C-rich case. That 
blending CN line seems to be real, as the shape of the two lines 
also differs in the observations, but the strength of the CN line 
is wrong in the synthetic spectrum (independent of changes in 
the global parameters). Thus this line could not be used and the 
determination of the isotopic ratio could be based only on one 
line with some additional indications from one blended feature 
in the C-rich case. 

3.5. Evolutionary models 

The theoretical models we use in the discussion of the observa- 
tional results have been obtained by means of the FRANEC stel- 



Table 1. Radial velocities of NGC 1846 AGB stars. Concerning 
the membership of LE17 please see the discussion in the text. 
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lar evolutionary code (Chieffi et al. |1998l l. A detailed description 
of the input physics can be found in Straniero et al. (2006) and 
Cristallo et al. (2007) . We calculated evolutionary sequences for 
various initial masses and compositions. In particular, models 
labeled as scaled solar, were obtained by scaling the solar abun- 
dances derived by Lodders(2003 ) of a factor Z/Zq, where Z is 
the assumed metallicity of the model. Then, a-enhanced models 
were obtained by increasing the initial abundance of O, Ne, Mg, 
Si, S, Ca, and Ti by a common factor with respect to the corre- 
sponding scaled solar values (see Piersanti et al. i2007l for further 
details). 

4. Results 

4.1. Cluster membership 

From the K-band spectra, we measured radial velocities of the 
target stars from individual lines of CO. Results are presented 
in Table[T] For LE8, Olszewski et al. ( I1991I I give a velocity 
of 248+5 km s"', in excellent agreement with our value. The 
mean velocity of all investigated AGB stars in NGC 1846 is 
248+5.5 km s"', which is in good agreement with the published 
values for the cluster. The only star with a larger deviation from 
the mean velocity is LE17. 

4.2. O-rich stars 

The C/O ratio and the '^C/'^C ratio were determined for all 
six O-rich stars in our sample. The C/O ratios between 0.2 
and 0.65 were found with an uncertainty between +0.05 and 
+0. 1, and '^C/'^C isotopic ratios between 12 and 60 were found. 
Temperatures of the best fit model spectra agree closely with the 
Teff values derived from near infrared photometry (Lebzelter & 
Wood 2007). The result for each star is listed in Table|2]with the 
basic parameters of the stellar models used. 

The C/O ratio obviously varies significantly from star to star, 
indicating different amounts of enriched material dredged up to 
the surface. It would be expected that an increased C/O ratio due 
to the dredge up of '^C would be accompanied by an increased 
isotopic abundance ratio. Indeed we find a nice correlation be- 
tween C/O and '^C/'^C for the O-rich stars as can be seen from 
Table|2] The values for C/O and '^C/'^C ratio found are in good 
agreement with findings from Smith & Lambert d 19901 1 for field 
stars. The minimum C/O ratio in our sample of 0.2 agrees closely 
with expectations from first dredge-up. 

As described above, one blend including a line of HF (1-0 
R7) is covered by our K band observations. After limiting the 
various parameters affecting the shape and depth of this blend, 
namely mainly the temperature and the carbon isotopic ratio, 
from other parts of the stellar spectra we fitted the remaining 
difference between the observed and the calculated blend profile 
by a change in the F abundance. Due to the rather large uncer- 
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Table 2. Resulting abundances and model fit parameters for the sample stars. The last two columns give the periods and luminosities 
taken from Lebzelter & Wood ( [2007] l. 
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Fig. 2. Part of the K band spectrum of LE8. The location of the 
HF line is marked by an arrow. Beside the spectrum (line + sym- 
bol) three models are plotted with a variable F abundance of -0. 1 
dex (dotted), -0.5 dex (solid), and -0.7 dex (dash-dot) relative 
to a solar value (-7.44) and scaled for the average cluster metal- 
licity, respectively. 



tainties of the C-star fitting we decided to derive the F abundance 
only for the O-rich stars. 

Figure |2] shows as an example part of the spectrum of LE8 in- 
cluding the HF line (marked by an arrow). Overplotted are three 
models with different F abundances. Obviously, Fluorine is un- 
derabundant in LE8 relative to the other metals. It turns out that 
the F abundance has to be changed from star to star to allow for 
a reasonable fit of the observed blend. Uncertainty in determin- 
ing the F abundance comes from both the temperature and the 
'^C/'^C ratio (via a '^CO line in the blend). We made an esti- 
mate of the uncertainty by varying these two parameters around 
the values determined from other parts of the spectra and then 
summing in quadrature the respective uncertainties. Taking into 
account this uncertainty and the fact that we have only one HF 
line to derive the abundance the results of course have to be taken 
with some caution. Confirmation from other HF lines is needed 
to limit the possible effect of unidentified lines on the profile of 
the blend and the derived abundances. 

4.3. C-rich stars 

The C/O ratio was derived for two C-rich stars, LEI 1 and LE2. 
For the other C-rich stars, either no H band spectra were avail- 
able (as a result of limited observing time) or no acceptable fit of 
the spectrum could be achieved. For the latter case, the starting 



influence of stellar variability on the stellar spectrum may play 
some role besides the above-mentioned problems. 

The derived T^ff and C/O for the C-rich star LEI 1 come with 
a much larger uncertainty. From the strength of the CO 3-0 band 
head, we exclude a C/O of less than 1 .4. A macroturbulent veloc- 
ity of 10 km s ' is needed to both fit the width of the features and 
their strength. The temperature can be constrained even less. If 
we use the value from near infrared photometry, we get a T^s of 
2950 K. As noted by Lebzelter & Wood (2007), this temperature 
is likely to be too low as pulsation properties instead suggest a 
value close to 3150K. We produced model spectra around both 
temperatures. A value of Teff =2950 K seems to give the better fit, 
but the sensitivity of the synthetic spectrum on this parameter is 
not very high. 

Concerning an upper limit of the C/O ratio an overall fit of 
the spectrum of similar quality can be reached for various com- 
binations of C/O and Teif with a higher temperature requiring a 
higher C/O ratio. As a rule of thumb an increase in T^ff by 100 K 
can be compensated by an increase of C/O by 0. 1 in the param- 
eter region under consideration. Taking the uncertainties of the 
line lists into account we cannot exclude a C/O ratio as high as 
2.0, but an upper limit of 1 .9 seems reasonable due to constraints 
set by the range of possible temperatures defined by the K band 
spectrum, the near infrared color and the pulsational behavior, 
respectively. In Fig. [3] we show a fit with C/0=1.7 and T^ff of 
2950 K. 

With this value we derived a '^C/'-'C ratio of about 60 from 
the K band spectrum. The derived isotopic ratio depends on the 
temperature, the C/O ratio, and the selection of the point in the 
synthetic spectrum to which the observed one is scaled. The C/O 
ratio has the least influence and, within the constraints set by the 
H band spectrum, its uncertainty can be neglected. The tempera- 
ture is an important parameter for Teff above 2900 K. Below that 
value, uncertainties in this parameter affect the resulting '^C/'^C 
ratio less and less. The selection of the 'pseudocontinuum' is of 
high importance in this low temperature regime. From a com- 
parison between the observed spectra and synthetic spectra by 
varying the mentioned parameters, we derive a possible range 
for the isotopic ratio of ±10 around the mean value. 

Additionally we were able to produce a reasonable fit for 
the spectra of LE2 when using a temperature between 2600 and 
2850 K, the H band spectrum indicates a C/O ratio of 1.9+0.2. 
The lower temperature nicely corresponds to the one derived 
from the J - K value (Lebzelter & Wood l2007b . The resulting 
'^C/'^C ratio is close to 65±15. 

While we were not able to make a good fit of the H and K 
band spectra of LEI, we made a qualitative comparison of LEI 
and LE2. Both wavelength ranges indicate that LEI has a higher 
C/O ratio than LE2. Unfortunately, a quantitative estimate was 
not possible. 
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Fig. 3. Spectrum of the carbon star LEll compared to a model 
fit with C/0=1.7 and reff=2950K. See text for details. 



5. Discussion 

5.1. Dredge up along the AGB 

In Fig.lHthe location of the O-rich stars of our sample in a color- 
magnitude diagram are shown. The data points are labeled with 
their respective C/O ratio. Color and m/,,,/ values used here are 
taken from Lebzelter & Wood (2007 ). As shown in that paper, 
color variations of these stars are very small and can be neglected 
in the present discussion. The mhoi values are derived from the 
near infrared magnitudes using the transforms in Houdashelt 
et al. (2000a 2000b). For details we refer to Lebzelter & Wood 
( |2007j . The sample stars cover a significant range in brightness 
along the AGB, thus allowing us to look for changes in the C/O 
and '^C/'-'C ratios with increasing luminosity. We can exclude 
that these stars belong to the RGB. Indeed, the RGB tip observed 
by Lebzelter & Wood ( 2007 ) is definitely fainter than the faintest 
star in the sample, LE17. This occurrence is confirmed by our 
evolutionary models. Actually, LE17 is only slightly brighter 
than the RGB tip of the 1.5 Mq model, but its '^c/'^c ratio (30) 
is definitely greater than expected for the brightest RGB stars 
(see e.g. Sneden fTWT] Gratton et al 2000). LEO, the brightest 
star in our sample, is also the brightest O-rich AGB star known 
in this cluster. However, the resulting picture is far from clear 
Indeed we find the star with the highest C/O value at the top. But 
the two stars with a C/O ratio close to 0.44 (LE16 and LE17) do 
not seem to fit into a steady increase with luminosity at all. They 
also seem to be offset from the other stars, which seem to form 
some sequence in the CMD. 

We note here that the luminosity at which an AGB star is ob- 
served does not necessarily have to be representative of its evo- 
lutionary state along the AGB. During thermal pulses the bright- 
ness can change significantly both to lower and to higher values. 
It was suggested in Lebzelter & Wood (2007 ) that at least the 
atypically low luminosity of two of the carbon stars in this clus- 
ter can be explained by a deviation from the mean interpulse 
luminosity during a thermal pulse. 

LE17 seems to be an outlier in various aspects. Assuming 
that the star is a cluster member, its relatively low luminosity, 
coupled to its relatively large J-K, could be explained if this 
star is undergoing a post-flash dip phase (Iben & Renzini [19831 1. 
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Fig. 4. Color-magnitude diagram of the O-rich stars in our sam- 
ple. The C/O ratio of each data point is given. 



The expansion powered by a thermal pulse indeed causes a sig- 
nificant decrease in the luminosity and in the effective temper- 
ature. From our models, we derive that variations up to 1 mag 
and 200 K (in Mboi and Tetf, respectively) have to be expected. 
Alternative scenarios for explaining the characteristics of LE17 
could be either an unresolved binary or a long secondary pe- 
riod similar to the ones already found in several AGB stars (e.g. 
Wood et al. 120041 ). The observed offset in velocity of LE17 from 
the other AGB stars in our sample could be consistent with both 
of these scenarios. The star's obvious enrichment in '^C, how- 
ever, could be best explained by the post-flash dip scenario. 

An alternative would be that LEI 7 is not a cluster member 
but a background star. Indeed this solution cannot be ruled out. 
Frogel et al. (1990) list the star as a cluster member based on 
its projected distance to the cluster center Lebzelter & Wood 
(120071 ) used the location of the star in the color-magnitude di- 
agram as an argument in favor of the cluster membership. On 
the other hand, the radial velocity of the star is offset from the 
other AGB stars measured by almost 20 km s ' . In their work on 
the velocity dispersion in LMC clusters Grocholski et al. (120061 ) 
chose a range in velocity of +10kms"' around the mean clus- 
ter velocity to separate members from non-members. According 
to this rule LEI 7 would not be a member Based on the present 
data, it is difficult to distinguish between the two possibilities. 

Concerning the other "anomalous" star, LE16, one could 
speculate that the observed AGB stars may belong to two se- 
quences, one including the stars with C/O ratios between 0.2 and 
0.3, and another one, shifted to the blue, with the two remaining 
stars. We come back to this point below. 

As mentioned above there is a qualitative agreement in the 
derived abundances between the increase in C/O and '^C/'^C. 
Thus, a wrong abundance for the two deviating stars is not very 
likely. Note that all the stars with C/O>0.2 and '^C/'^OIS but 
LEI 7 are brighter than the minimum luminosity for the occur- 
rence of the third dredge up, as predicted by our models. Taking 
the mean luminosity of the interpulse period just after the occur- 
rence of the first TDU episode, we find that all stars with mhoi 
brighter than 14.1 mag should be '^C enhanced. Once again, the 
anomaly of LE17 could be solved if the star was in the post-flash 
dip. 

We illustrate flie relation between C/O and '^C/'^C in Fig.|5] 
While there is some scatter, a relation between these two quan- 
tities can clearly be seen. In the same figure, we also plotted the 
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Fig. 5. C/0 versus '^C/'^C for the O-rich stars in our sample. For 
comparison predicted values from our evolutionary models are 
shown. The lines connect individual interpulse values for four 
sets of model parameters. See text for details. 



expected relations from evolutionary models. The solid line rep- 
resents our reference model: M=1.9 M© (1.8 at the beginning 
of the TP- AGB phase), Z=0.006, [a/Fe]=0.2, Y=0.27. The mass 
has been chosen according to the pulsational mass derived by 
Lebzelter & Wood iQQGJj, while the moderate enhancement of 
the a elements corresponds to the value generally claimed for the 
Large Magellanic Cloud (see e.g. Hill et al. 2000). The C/O ratio 
of the two less evolved stars in our sample, H39 and LE9, sup- 
ports such an assumption. Indeed, the expected C/O ratio after 
the first dredge up is 0.35, for a scaled solar composition, and 0.2 
when [OfPe]-+0.2. In this case, the total metallicity (Z=0.006) 
would imply [Fe/H]=-0.5, in agreement with the latest spectro- 
scopic determination for NGC 1846 (Grocholski et al. 120061 ). 
To evaluate the effects of possible variations of these model pa- 
rameters, we also report the relations obtained by changing the 
mass (1.5 Mq, dashed fine), the total metallicity (Z=0.003, dot- 
ted line), and the initial abundance ratios (scaled solar instead of 
a-enhanced, dot-dashed line). 

It can be seen that all data points are found very close to 
the theoretical relation for the reference model, even if LEI 7 
presents a lower carbon isotopic ratio (or too high C/Ofl 
Changes in initial mass and metallicity have little effect, so that 
the uncertainties of the present spectroscopic analysis do not al- 
low us to distinguish among the different relations reported in 
Fig. 5. The model with a solar-scaled composition can instead 
be safely excluded. The boxes along the relations mark the inter- 
pulse values after each dredge up event. Close to these values we 

^ Note that the abundances have been determined with a model at- 
mosphere without O enhancement. To see if this has any effect on the 
derived ratios, we calculated a test model with O enhanced by +0.1 
and +0.2 dex, respectively. It turned out that this enhancement only has 
a weak effect on the spectra which can be compensated for by small 
changes in the temperature. The C/0 and the '-C/'^C ratios remain al- 
most unchanged. 



should have the highest probability of finding a star. In general, 
the agreement is remarkably good, with one exception (LEO). 

Concerning the abundance of the light element F, we ob- 
serve, as shown in Fig.|6] a change from an F underabundance 
to an F overabundance with increasing C/0 (relative to the solar 
value scaled to the iron abundance of the cluster). Accordingly, 
there is also a similar trend with the luminosity, as found for the 
C/O ratio. Measurements of HF lines are a critical source for 
obtaining the F abundance in stars. Various hypotheses on the 
origin of this element have been discussed in the literature (see 
Schuler et al. 2007 for a recent overview). In particular, from 
a correlation of the carbon abundance with the F abundance, 
Jorissen et al. ( 119921 ) showed that nucleosynthesis in AGB stars 
is one likely path for the production of this element. Schuler et 
al. (2007) come to a similar result from F abundance measure- 
ments of a carbon-enhanced, metal-poor star A positive correla- 
tion between F and C/O has also been acknowledged in a plan- 
etary nebula (Zhang & Liu "SOOS). Our data of NGC 1846 AGB 
stars clearly support the production and dredge up of F during 
the AGB evolution. The two lines in the figure represents the 
predictions of the reference model. The only difference between 
the two lines is in the initial fluorine abundance. We have as- 
sumed [F/Fe]=0 and [F/Fe]=-0.71 at the beginning of the evolu- 
tionary sequence, for the solid and the dashed lines, respectively. 
The lower value corresponds to the one of the less evolved AGB 
star in our sample (H39). Note that Cunha et al. (2003) report 
similar fluorine underabundances in LMC giants. Our models 
predict a negligible variation in the surface F abundance after 
the first dredge up. It should be noted that only an "extreme" 
extra mixing during the RGB may induce a sizeable F deple- 
tion (see e.g. Denissenkov et al. .2006 ). In contrast, during the 
AGB fluorine is produced by a captures on '^N in the convec- 
tive zone generated by a thermal pulse and, later on, dredged 
up. Part of the '^N seeds are in the ashes of the CNO burning 
(about 50%), while the remaining part is synthesized in the '^C 
pocket. In both cases, '^N is released by the '^0(/:', q')'^N reac- 
tion (for more details see Cristallo et al. l2006l ). However, even 
if models account for a substantial AGB production, we confirm 
the previous claim by Jorissen et al. ( 1992) that the theoretical 
predictions underestimate the fluorine overabundance observed 
in evolved AGB stars. A possible solution may be a more effi- 
cient production of '^N in the He-rich inter-shell zone. Perhaps, 
an additional and different source of fluorine should be sought. 
We have to note, of course, that, as the fluorine abundance in the 
NGC 1846 stars has been derived from only one line blend, a 
systematic error cannot be excluded. 

5.2. Dredge up and pulsation properties 

The detailed variabiUty study on the A GB stars in NGC 1846 
undertaken in Lebzelter & Wood ( 120071 ) allows us to investigate 
the observed star-to-star differences in the C/O ratio in light of 
the pulsational properties of these stars. In Fig.|7]we show the lo- 
cation of our sample stars in the P-L-diagram of the cluster. Data 
are taken from Lebzelter & Wood. LE9 could not be plotted due 
to the lack of any detected periodicity in its small-amplitude light 
variability. The other stars can be easily attributed to a pulsation 
mode: LEO and LE16 are first overtone pulsators, while LE8 
and H39 pulsate in the second overtone. Long secondary peri- 
ods beyond the fundamental mode period have been found in the 
light curves of several of these stars, but no shorter secondary pe- 
riods. As already pointed out in Lebzelter & Wood (120071 ). LE17 
cannot be attributed unequivocally to any pulsation sequence. 
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Fig. 8. Location of the O-rich stars in a color-magnitude diagram 
with the pulsation mode of each object indicated. Open boxes 
denote first overtone pulsators, while filled boxes mark stars with 
second overtone pulsation. The star LE17 has been excluded. 



Fig. 6. Fluorine abundances of our O-rich stars sample as a func- 
tion of the C/O ratio. The solid line corresponds to the reference 
theoretical model, while the dashed line represents a model for 
which the initial [F/Fe] is -0.7 1 instead of 0. 
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Fig. 7. The logCL/L©) - log P diagram for NGC 1846 adapted 
from Lebzelter & Wood ( I2007I I. Pulsation models for fundamen- 
tal mode and the first three overtone modes are shown together 
with five O-rich stars of our sample. 



Stars located on the first overtone sequence show a higher 
C/O ratio than objects on the second overtone sequence (Fig.|2l). 
As can be seen this cannot be attributed to a simple difference in 
luminosity as the two groups clearly overlap. This is further il- 
lustrated in Fig.[8]where we marked the sample stars in the color- 
magnitude diagram according to their pulsation mode. For LE9 
(C/O=0.2), no periodic change was found by Lebzelter & Wood 
( 1200 7). Assuming that it would have been more likely to find 
a first overtone pulsation than a shorter period higher overtone 
pulsation, we attributed this star to the group of higher overtone 
pulsators. These findings suggest that in addition to a relation 



between the C/O ratio and the luminosity, seems to exist some 
dependency of this ratio on the pulsation mode. At this point it 
is not clear how to explain the existence of such a relation. As 
mentioned above, all sample stars are on the AGB, we can ex- 
clude an RGB contamination. From the difference in the C/O 
ratio, we can assume that the sample stars are in at least three 
different thermal pulse cycles. Both from our models and from 
the literature (e.g. Vassiliadis & Wood 1993 ), we would expect 
that the stars become cooler from cycle to cycle. As the two C 
enhanced stars LEI 6 and LEO would form a sequence at higher 
temperatures than the other stars, it can be excluded that a differ- 
ence in the thermal pulse cycle is responsible for the observed 
differences in pulsation mode and temperature. 

A difference in pulsation mode between two stars of the same 
luminosity requires a difference in some other global character- 
istic of the star. This could be a difference in mass. Mackey & 
Broby Nielsen ( I2007I I report indications for two stellar popula- 
tions in NGC 1846 separated in age by about 300 Myr. At a 
given time we may then find AGB stars of different masses at the 
same time on the AGB. Indeed such an age or mass difference 
could easily account for the observed color difference (Fig.lSj. 
This would also be compatible with the observed C/O ratios as 
the more massive stars may show a more efficient dredge up. 
Furthermore, the starting abundances of the two populations may 
be different, which could enhance the observed difference in C/O 
and i^c/i^c. 

It has to be stressed that, even in the case of a single stellar 
population, a scatter in the abundances of C and O on the main 
sequence and after the first dredge up may be present. Such an 
abundance scatter has been found in various clusters. 

5.3. The S-type stars 

An interesting case is LEO. This star has been classified as 
spectral type S3/3 by Lloyd Evans (1983i ). Indeed the ZrO band 
head at 6345 A is clearly visible in the spectrum presented by 
him (see his Fig. 1). According to its spectral classification, TiO- 
bands and ZrO bands should be of similar strength (Keenan & 
Boeshaar 119801 ). YO should occur, but is not covered by the 
spectra shown by Lloyd Evans. As expected for a star with S- 



T. Lebzelter et al.: AGB stars of the intermediate-age LMC cluster NGC 1846 



9 



T ' 1 ' T 




wavelength [A] 

Fig. 9. H band spectrum of LEI 3 compared with models of 
3600 K, logg=0, C/O=0.6 (solid line) and 0.7 (dotted line), re- 
spectively. 



type spectral characteristics, the C/O ratio is clearly enhanced in 
LEO with a value between 0.6 and 0.7. In Fig.|9]we compare 
the observed spectrum with our model for both cases (C/O=0.6 
and C/O=0.7). The model obviously fits the observations quite 
well. As mentioned above, this model was calculated with an in- 
creased microturbulent velocity to allow for a reasonable fit of 
the observed line profiles. A change in the C/O ratio to 0.7 gives 
a similar fit, but a further increase in the C/O ratio gets problem- 
atic as the CN lines, especially the one at 15567 A, become too 
strong. Compensating for this with a higher temperature would 
give a bad fit of the OH line close to 15570 A and too strong a 
CO band head. Similarly, changes in log g cannot compensate 
for a C/O value close to one. LEO is thus an S-type star with 
C/O clearly below 1 . 

Lloyd Evans ( 1983) also classified a second star in this clus- 
ter as of spectral type S3/1, namely LE8. LE8 has an even lower 
C/O ratio of only about 0.3. Direct comparison with the spec- 
tra of LE9 and H39 does indicate that the star has an increased 
C/O ratio, but it is far below the expectations for an S-star. 
Unfortunately, Lloyd Evans (I 1983I I presented no visual spectrum 
of this star. 

A C/O ratio of only 0.65 for an S-type star is in agreement 
with the oxygen overabundance in the Magellanic Cloud dis- 
cussed above. If this the case for NGC 1846, and the comparison 
of the other stars with the expectations from the model suggests 
this, then we would require more dredge up events to produce a 
C/O ratio close to one than in the case of a solar abundance pat- 
tern. In other words, the abundance of s-process enriched mate- 
rial on the surface may be high enough to see e.g. ZrO bands, but 
the amount of C dredged up does not yet equal the amount of O 
due to the high starting value of the O abundance. Comparing the 
C/O ratio of LEO and LE8, we would expect to see less ZrO in 
LE8, which indeed agrees with the spectral classes from Lloyd 
Evans (1983). 

Piccirillo d 19801 1 pointed out that, only for stars with 
T<3000 K, would we expect that the occuiTence of ZrO bands 
correlates with a C/O close to one. For hotter stars, an increased 
abundance of Zr is the main reason for a visibility of ZrO in the 
spectrum. As the stars in our sample all have temperatures above 
3000 K, the occurrence of ZrO bands is not surprising. We think 



that a combination of this effect and an enhancement of O in the 
pre- AGB composition of the star explains the observed C/O ra- 
tio. In this context we also note that Smith & Lambert (199Q| 
report a C/O for the two intrinsic (i.e. AGB, see Jorissen et al. 
1993) S-stars NQ Pup and V679 Oph of 0.29 and 0.75, respec- 
tively. Thus the C/O values we found for the two S-type stars in 
NGC 1846 look realistic. 



5.4. The C-stars 

The study of the C-stars in this cluster was hampered by these 
difficulties deriving basic parameters and abundances. Still a few 
points can be noticed. The C/O ratio we find is rather high. There 
is quite a gap in our sample between the highest C/O ratio of an 
O-rich star (LEO, 0.65) and the lowest value of a C-rich star 
(LEll, 1.7). There are still two candidates in the cluster, LE12 
and LW2 (see Lebzelter & Wood 2007), which may fill this gap, 
but unfortunately we have no data for these stars. The value is 
also high compared to the values found by Lambert et al. ( 119861 
see also Smith & Lambert fl 9901 1 for field stars, where the maxi- 
mum C/O is 1.6. Such an occurrence is not surprising, however, 
since the oxygen content is lower in NGC 1846, even consider- 
ing [O/Fe]=0.2, when compared to field stars. From the abun- 
dance of HCN and C2H2, Matsuura et al. (2006) find that the 
C/O ratio in their sample of LMC carbon stars should be higher 
than 1.7, which would agree with our findings. 

Of particular interest is the carbon isotopic ratio derived for 
the 2 carbon stars in our sample. Although the C/O ratio in- 
creases, the '^C/^-'C remains close to the highest value found 
for the O-rich stars in this cluster. Basing on the only 2 C stars 
in the sample, we can tentatively conclude that some saturation 
level exists around 60 for the carbon isotopic ratio. In contrast 
to this figure, the evolutionary models predict a continuous in- 
crease in the '^C/'-'C ratio, as due to the dredge up of primary 
'^C. Indeed, for C/0> 1 , our reference model predicts values def- 
initely higher than 100. Our finding confirms what has akeady 
been observed in galactic C-stars (Smith & Lambert 1990. Abia 
et al. 120011) . which shows a clustering between 50 and 70. 

The most obvious explanation for this observational pattern 
is the occurrence of a mixing process able to bridge the radiative 
gap between the cool bottom of the convective envelope and the 
hot H-burning zone. A similar process seems to explain some 
abundance anomalies in red giants stars (see, e.g., Charbonnel 
1 19951 ). It is usually referred to as extra or deep mixing. More 
recently, NoUett et al. (2003) have also proposed a deep circu- 
lation in low mass AGB stars as a solution to the puzzling iso- 
topic anomalies found in some pre-solar grains. They named this 
circulation the cool bottom process. In spite of all the observa- 
tional evidence, a common consensus on the physical mecha- 
nism driving this mixing has not been established yet (e.g. ro- 
tational induced instabilities, magnetically induced circulation, 
gravity wave, or thermohaline mixing - see e.g. Busso 120071 for 
a recent review). 

The scenario emerging from the evolutionary sequence of 
'^C/'^C versus C/O in NGC 1846 is that of a moderate deep mix- 
ing, affecting the '^^C surface abundance in the late part of the 
AGB only, when stars become C-rich. Indeed, while the abun- 
dance derived from O-rich stars in NGC 1846 are all compatible 
with the prediction of models with no extra mixing, the '^C/'^C 
ratios we find for the two C stars indicate that their envelope ma- 
terial could have been exposed to a temperature on the order of 
30-40 ■ lO*" K. At that temperature and on an AGB timescale, '^C 
is partially converted into '^'C, but only marginally into '"^N. 
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Fig. 10. C/O versus '^C/'^C for all the sampled stars. Predicted 
values from models with and without extra-mixing are shown. 
Solid line refers to our reference model, while the dotted 
and dashed ones refer to models with T/,m=35-10'' K and 
T/,m=40 -10^ K, respectively (see text for details). 



From the models, we find that, during the AGB the bottom 
of the convective envelope, whose temperature never exceeds 
4 ■ 10'' K, moves closer to the H-burning zone. Just before the 
first thermal pulse, a region as large as 2.2 ■ 10""' M© separates 
the innermost convectively unstable layer from the shell where 
the temperature is 40 ■ 10^ K, but before the 10th TP, this region 
is reduced to 4 ■ 10"** M©. Such an occurrence suggests that the 
effect of a possible extra mixing should be stronger in the late 
part of the AGB than at the beginning. 

To simulate such a scenario, we calculated an additional 
model with the same input parameters as the reference model, 
but artificially including an extra-mixing. We assume that such 
an extra-mixing is confined within a region as large as 10"^ M© 
below the convective envelope, but only if the temperature is 
lower than a maximum value (T/„„). The average convective ve- 
locity is fixed to 100 cm s"' 0. The most important parameter is, 
here, the value of the maximum temperature, as shown in Fig.fTOl 
where we reported the '^C/'^C versus C/O relations obtained for 
T/,„,=35 and 40 (10'' K). 

This phenomenological model is useful for inferring the ex- 
tension of the region affected by the supposed extra-mixing and 
its maximum temperature, but a more physical description of 
the mechanism supporting this phenomenon is needed to derive 
more quantitative results. Either way, we stress that the '^C/'^'C 
versus C/O sequences for intermediate age clusters are the best 
way to understand and to constrain the deep mixing in AGB 
stars. 



^ This is a slow mixing, when compared to the typical convective 
velocity in the envelope of an AGB star, which is of the order of 10^ 
cm/s. We have verified, however, that the C/O and '^C/'^C ratios are 
only marginally affected by the choice of this parameter. 



6. Conclusions 

We derived C/O ratios and '^C/'^C ratios for the first time for 
a sample of AGB stars above the dredge-up mass limit within 
one stellar cluster In this way we tested how the various third 
dredge-up events change the surface abundances in a homoge- 
neous sample of AGB stars. Oxygen-rich stars with clear indi- 
cations of a third dredge-up could be identified from enhanced 
values of these two quantities, which lie between 0.2 and 0.65 
for C/O and 12 to 60 for '^C/'^C. The cluster also contains a 
few C-stars, but only rather uncertain values for the C/O ratio 
could be derived. A more detailed analysis for two C-stars was 
possible leading to C/0=1.7 and 1.9 and i2C/'^C=60 and 65. 

As expected, the C/O ratio is correlated with the '^C/'^C ra- 
tio. We find good agreement with model predictions for both ra- 
tios assuming an oxygen overabundance of +0.2 dex in the clus- 
ter. The measured C/O ratios are in reasonable agreement with 
the predicted interpulse values. This is the first direct test of the 
models concerning dredge up before the stars become C-stars. 

A possible correlation between C/O ratio and pulsation mode 
was found. The difference in pulsation mode may also help 
to explain the observed scatter in C/O along the AGB of the 
cluster, although no final understanding of this relation could 
be achieved. The findings may be related to the existence of 
two populations within the cluster with a small difference in 
age. Further investigations of the populations in this cluster are 
needed. 

An increase in the fluorine abundance along the AGB has 
been found. This confirms nucleosynthesis in AGB stars as one 
of the origins of this element. 

Two stars in the cluster were classified as having spectral 
type S in the past. We find a C/O of 0.65 and 0.3 for these objects. 
While the C/O ratio is thus clearly less than 1, we argue that 
the comparably high temperature of the stars allows for S-type 
characteristics at this C/O. 

The '^C/'^C ratio found in the C-stars cannot be explained 
with standard models for the third dredge up. They indicate the 
need for a moderate extra-mixing, affecting the late part of the 
AGB evolution. 
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